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In this paper, we obtain higher dimensional topological black hole solutions of Einstein-A gravity 
in the presence of a class of nonlinear electrodynamics. First, we calculate the conserved and 
thermodynamic quantities of (n + l)-dimensional asymptotically flat solutions and show that they 
satisfy the first law of thermodynamics. Also, we investigate the stability of these solutions in the 
(grand) canonical ensemble. Second, we endow a global rotation to the static Ricci-flat solutions 
and calculate the conserved quantities of solutions by using the counterterm method. We obtain a 
Smarr-type formula for the mass as a function of the entropy, the angular momenta and the electric 
charge, and show that these quantities satisfy the first law of thermodynamics. Then, we perform a 
stability analysis of the rotating solutions both in the canonical and the grand canonical ensembles. 


I. INTRODUCTION 

Nonlinear field theories are of interest to different branches of mathematical physics because most physical systems 
are inherently nonlinear in the nature. The main reason to consider the nonlinear electrodynamics (NLED) comes 
from the fact that these theories are considerably richer than the Maxwell field and in special case they reduce to the 
linear Maxwell theory. Various limitations of the Maxwell theory, such as description of the self-interaction of virtual 
electrompositron pairs and the radiation propagation inside specific materials .3f], motivate one to consider 

NLED Q. Besides, NLED improves the basic concept of gravitational redshift and its dependency of any background 
magnetic field as compared to the well-established method introduced by standard general relativity. In addition, it 
was recently shown that NLED objects can remove both of the big bang and black hole singularities Moreover, 
from astrophysical point of view, one finds that the effects of NLED become indeed quite important in superstrongly 
magnetized compact objects, such as pulsars and particular neutron stars (also the so-called magnetars and strange 
quark magnetars) Q. 

About eighty years ago Born and Infeld introduced an interesting kind of NLED in order to remove the divergence 
of the self-energy of a point-like charge j7|. The first attempt to couple the NLED with gravity was made by Hoffmann 
(§]. After that the effects of Born-Infeld (BI) NLED coupled to the gravitational field have been studied for static 
black holes Q , rotating black objects 0 , wormholes [HI, [HI and superconductors 0. Also, BI NLED has acquired 
a new impetus, since it naturally arises in the low-energy limit of the open string theory [l4| . Recently, two different 
BI types of NLED have been introduced, which can also remove the divergence of the electric field near the origin. 
One of them is Soleng NLED which is logarithmic form ITbl and another one was proposed by Hendi with exponential 
form 0], The Soleng field, like BI theory, removes divergency of the electric field, while the theory proposed by 
Hendi does not. It is notable to mention that although the exponential form of NLED does not cancel the divergency 
of the electric field but its singularity is much weaker than that in the Maxwell theory. Black object solutions coupled 
to these two nonlinear fields have been studied in literature (for e.g., see 00). The Lagrangian of mentioned BI 
type nonlinear theories, for weak nonlinearity, tends to the following form 

L{F) = —T + afF 2 + O (a 2 ), (1) 

where T = F IJl/ F ,JL ' is the Maxwell invariant, in which F^ = is the electromagnetic field tensor and A^ 

is the gauge potential. In addition, a denotes nonlinearity parameter which is small and so the effects of nonlinearity 
should be considered as a perturbation (a is proportional to the inverse value of nonlinearity parameter in Bi-type 
theories). In this paper, we take into account the Eq. (JT]) as a NLED source and investigate the effects of nonlinearity 
on the properties of static and rotating black hole/brane solutions. 

Here, it is necessary to focus on the basic motivation of considering the Lagrangian ©• At first, we should note 
that, regardless of a constant parameter, most of NLED Lagrangians reduce to Eq. © for the weak nonlinearity. 
Eventually, it is worthwhile to mention that although various theories of NLED have been created with different 
primitive motivations, only for the weak nonlinearity (Eq. ©), they contain physical and experimental importances. 
As we know, using the Maxwell theory in various branches leads to near accurate or acceptable consequences. So, 
in transition from the Maxwell theory to NLED, the logical decision is to consider the effects of weak nonlinearity 
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variations, not strong effects. This means that, one can expect to obtain precise physical results with experimental 
agreements, provided one regards the nonlinearity as a correction to the Maxwell field. 

For the reasons mentioned above, there have been published some reasonable works by considering Eq. © as an 
effective Lagrangian of electrodynamics 1, 2, 4, 18H22|. Heisenberg and Euler have shown that quantum corrections 
lead to nonlinear properties of vacuum .1, [2| j_4J, [l8|. Also, as we mentioned before, it was proved that in the low- 
energy limit of heterotic string theory, a quartic correction of the Maxwell field strength tensor appears (Til - IEil . So it 
is natural to consider Eq. © as an effective and suitable Lagrangian of electrodynamics instead of the Maxwell one. 

The outline of our paper is as follows. In the next section, we consider the (n + l)-dimensional topological static 
black hole solutions of Einstein gravity in presence of the mentioned NLED and investigate their properties. In Sec. 
E© we calculate the conserved and thermodynamic quantities of asymptotically flat black holes, check the first law of 
thermodynamics and investigate the stability of the solutions in both canonical and grand canonical ensembles. Sec. 
II VI is devoted to introducing the rotating solutions with flat horizon and computing the conserved and thermodynamic 
quantities of the solutions. We also check the first law of thermodynamics and perform the stability analysis of the 
solutions both in the canonical and the grand canonical ensembles for the rotating solutions. We finish our paper 
with some concluding remarks. 


II. STATIC TOPOLOGICAL BLACK HOLE SOLUTIONS 

The (n+ 1)-dimensional action of Einstein gravity with negative cosmological constant and in presence of nonlinear 
electrodynamics is 


t g = T7 f d n+1 Xy/^lR-2A + L(F)}~±- [ d n xV^&h), ( 2 ) 

lOTT J M 8tT J dM 

where R is the scalar curvature, A is the cosmological constant which is equal to —n (n — 1) /2 1 2 for asymptotically 
adS solutions. In this action, L(J-) is the Lagrangian of nonlinear electrodynamics presented in Eq. © and the second 
integral is the Gibbons-Hawking surface term which is chosen such that the variational principle will be well defined 
[23 |. In the second integral, 7 and 0 are, respectively, the trace of induced metric, 7 ^, and the extrinsic curvature 
0,-j on the boundary dM. Variation of the action © with respect to the metric tensor g and the Faraday tensor 
F /JV . leads to 


G, 


nv + A g^ = -g^LiR) - 2L J rF fl \F l 




( 3 ) 


d» (V^gLjrF^) = 0 , ( 4 ) 

where G ^ is the Einstein tensor and Lj: = dL(T)/dT. 

Here, we want to obtain the (n + l)-dimensional topological static black hole solutions. We take into account the 
metric of (n + l)-dimensional spacetime with the following form 

dr 2 

ds 2 = -f{ r )dt 2 + ^- + r 2 dnl_ 1 , ( 5 ) 

9{r) 

where f{r) and g(r) are two arbitrary functions of radial coordinate which should be determined and 


n— 1i—l 

dQ 1+ © II sin 2 djdO 2 
*=2 j =1 

( n— 1 2—1 ^ 

dd2 + © II sin 2 djdO 2 

i=3 j =2 , 

n— 1 

E dti 
1=1 


k = 1 
k = — 1 , 
k = 0 


( 6 ) 


represents the line element of an (n — l)-dimensional hypersurface with constant curvature (n — 1 ) (n — 2 ) k and 
volume V n -\. We should note that the constant k characterizes the hypersurface and indicates that the boundary 
of t = constant and r = constant can be a positive (elliptic), zero (flat) or negative (hyperbolic) constant curvature 
hypersurface. 
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FIG. 1: elliptic horizon solutions: /(r) versus r for k — 1, n — 4, q = 1, A = —1, a = 0.005, and m = 1.1 (solid line), 
m = 1.2 (bold line) and m = 1.3 (dashed line); ’’dotted line is f(r) for the Maxwell case (a = 0) with m = 1.3 ”. 


Using Eq. J4]) with the following radial gauge potential ansatz 


A^ = h (r) <5°, 


( 7 ) 


we obtain the following differential equation 


E\r) 


(n — 1 )E{r) 
r 



1 )E(r) 


+ 3 E'{r) 


a = 0, 


( 8 ) 


where E(r) = F tr = h'(r ) is the nonzero component of electromagnetic field and prime denotes the derivative with 
respect to r. Solving Eq. ©, one obtains 


E(r) 


qa 2 r 2n+1 


{ 1 + y / 1 + i §0 


2/3 


— ar 


2 n 


Gar 71 


qa 2 r 2n+1 


( x + \/ 1 + iS') 


1/3 


(9) 


where q is an integration constant which is related to the electric charge of the black hole. Now, we use the series 
expansion of E{r) for small values of a, and keep the first two terms to obtain 


E(r) = 



or correspondingly 


h(r) = 


q 

(n — 2) r n ~ 2 


4 q 3 

(3 n — 4) r 3n ~ 4 


+ O (a 2 ) . 


( 10 ) 


( 11 ) 


It is easy to see that the second term in Eqs. m and m comes from the nonlinear correction and for vanishing 
a, one can reproduce the results of the Maxwell theory. Since we want to investigate the nonlinearity parameter as a 
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FIG. 2: flat horizon solutions: /(r) versus r for k = 0, n = 4, q = 1, A = —1, a = 0.02, and m = 0.4 (solid line), m = 0.5 
(bold line) and m = 0.6 (dashed line) ; ’’dotted line is f(r) for the Maxwell case (a = 0) with m = 0.6 ”, 


(perturbative) correction, hereafter, we take into account the first correction term of nonlinearity parameter, a, and 
ignore a 2 and higher power of nonlinearity parameter terms. To obtain the metric functions / (r) and g(r), one may 
use the nonzero components of Eq. ©• Straightforward calculations show that the nonzero components of Eq. © 
(up to the first order of a) can be written as 


ei = (n - 1 )rf(r) + (n - 1 )(n - 2 ) [f(r) - k] + 2Ar 2 + + O (a 2 ) = 0 , 

e 2 = r 2 f"(r ) + 2 (n - 2)rf'(r) + (n- 2)(n - 3) [/(r) - k] + 2 Ar 2 - + O (a 2 ) = 0. 

After some calculations, one can show that the solutions of Eqs. m and m can be written as 


( 12 ) 

(13) 


m 2A r 2 2 q 2 

r n -2 n ( n _ x) (n— 1) (n — 2 ) r 2n ~ 4 

_ 4g 4 _ 4- O ( 2 ) 

[2 (n - 2) (n + 2) + (n - 3) (n - 4)] r 4 "~ 6 ^ ’ ’ 


(14) 


ff(r) = C'/W, (15) 

where m is an integration constant which is related to the mass of the black hole and the last term in Eq. m 
indicates the effect of nonlinearity. Hereafter, we set the constant C = 1 without loss of generality. It is notable to 
mention that, for a = 0, this metric function reduces to Reissner-Nordstrom solution, as it should. The asymptotical 
behavior of the solution (fl4l) is adS or dS provided A < 0 or A > 0 and the case of asymptotically flat solutions is 
permitted for A = 0 and k = 1. 

Now we look for the singularities of the solutions. One can show that the metric © with the metric function (fldll 
has an essential singularity at r = 0 by calculating the Kretschmann scalar, as 


D rjjUZ'Att 

-tt/ii/A/c-nT 


8A 2 16Ag 2 (n 2 - l) (n 2 - 2) m 2 

n(n + 1) n{n + 1 )r 2n + r 2 n +2 

8 (2 n 2 — n — 2) mq 2 ( 96Aa 8 (8 n 2 — 12 n + 3) 

r 3n+i + ^ n + i+ n — 1 

32 (n — 1) (4n 2 — n — 2) maq 4 32 (32 n 2 — 32n + 5) aq 6 , 2 

(3n—l)r 5n+1 n(3n—l)r 6n + O [a 



( 16 ) 
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FIG. 5: Carter-Penrose diagram for the asymptotically flat (left figure) and the asymptotically adS (right figure) black holes 
when the metric function has one real positive root (ri) (the same as Schwarzschild black hole). 


r =0 



r =0 



FIG. 6: Carter-Penrose diagram for the asymptotically flat (left figure) and the asymptotically adS (right figure) black holes 
when the metric function has two real positive roots (ri and ?' 2 ) (second root is an extreme root). 


From Eq. © it is obvious that Kretschmann scalar diverges at r = 0 and, like the asymptotically adS solutions, 
it reduces to 8A 2 /n(n + 1) for r —> oo. 

Figs. |T]-[|]show that the singularity may be covered with horizon and, therefore, we can interpret the singularity as 
a black hole. In addition, these figures confirm that the nonlinearity parameter not only modify the electromagnetic 
part of solutions, but also the kind of horizons. For vanishing a, with suitable choice of parameters, metric function 
could acquire at most two horizons whereas, for this nonlinear theory (nonzero a), it is possible to find three horizons. 
Fig. |4] indicates that the nonlinearity parameter may considerably affect the existence, location and type of horizons. 
From Figs. |T]-[I]we depending the values of metric parameters with suitable a, the horizons of the black 

hole solutions may be extreme or not. 

In order to study the conformal structure of the solutions, one may use the conformal compactification method to 
plot the Carter-Penrose (conformal) diagram (see Figs. 07 ]). The Carter-Penrose diagrams and also the figures of 
the metric function (Figs. [T]-[3]) confirm that, the singularity is spacelike such as that of Schwarzschild black holes. 
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FIG. 7: Carter-Penrose diagram for the asymptotically flat (left figure) and the asymptotically adS (right figure) black holes 
when the metric function has three real positive roots (ri, ?'2 and 7 - 3 ). 


In other words, keeping the first order of nonlinearity parameter and ignoring the higher order of a, the timelike 
singularity of the Reissner-Nordstrom black holes (dotted line in Figs. [TJOI]) change to a spacelike singularity. Drawing 
the Carter-Penrose diagrams shows that the causal structure of the solutions are asymptotically well behaved. 

The temperature may be obtained through the use of regularity of the solutions at r = r+, yielding 


T+ = 


f(r+) 

47T 


1 


( (n — 1 ) (n — 2 ) k q 2 2 g 4 \ , 2 , 

(n — 1) l 2r+ + ri n " 3 rt n ~ 5 1 { ’ 


+ + / 

The electric potential d>, measured at infinity with respect to the horizon, is defined by {asm 

^ = |r—>-oo |r=r_|_? 

with the following explicit form 

q 4g 3 


4 > = 


( n — 2) 


„n —2 


(3n — 4) r\ 


3n—4 
+ 


oc T O . 


(17) 


(18) 


(19) 


The last term in the right hand side of Eqs. m and (flTH) indicates the nonlinearity effect of the mentioned NLED. 


III. THERMODYNAMICS OF ASYMPTOTICALLY FLAT BLACK HOLE (A = 0, k = 1) 

At first, we calculate the conserved and thermodynamic quantities of the black hole for A = 0 and k = 1. Second, 
we obtain a Smarr-type formula for the mass as a function of the entropy and the electric charge of the solutions and 
finally check the first law of thermodynamics. 
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The first quantity which we are going to calculate is the entropy of the black hole. More than thirty years ago, 
Bekenstein argued that the entropy of a black hole in Einstein gravity is a linear function of the area of its event 
horizon, which is the area law [HI]. Therefore, the entropy per unit volume V n -\ of the presented black hole is equal 
to one-quarter of the area of the horizon 


5 = 


( 20 ) 


In order to obtain the electric charge per unit volume V n -\ of the black hole, we use the flux of the electric field at 
infinity, yielding 

Q =b (21) 

which shows that, this kind of nonlinearity does not change the electric charge. The ADM (Arnowitt-Deser-Misner) 
mass of black hole can be obtained by using the behavior of the metric at large r. The mass per unit volume V n -i of 
the black hole is 


M = 


(n — 1 ) to 
167T 


( 22 ) 


where we can obtain m from f(r = r+) = 0. 

After calculating all of the conserved and thermodynamic quantities of the black hole solutions, we want to investi¬ 
gate the first law of thermodynamics. To do this, we obtain the total mass M as a function of the extensive quantities 
Q and 5. Using the expression for the entropy, the electric charge and the mass given in Eqs. mu, m and (1221) . 
and the fact that f(r = r + ) = 0, one can obtain a Smarr-type formula as 


M(S,Q) 


(n- 1) T r 
167r 


27 tQ 2 


64t r 3 Q 4 


(n - 2) T"- 2 (3n - 4) T 3i 


—a + 0(a 2 ) 


(23) 


where T = (4S') 1 ^” -1 ' ) . 

Now, we regard the parameters Q and S as a complete set of extensive parameters and define the intensive 
parameters conjugate to them. These quantities are the temperature and the electric potential 


dM\ _n{n-l){n-2)T n ~ 2 2i tQ 2 64t r 3 Q 4 . 2 , 

dS ) q 167m (n—ljS” 5 (n — 1) T n_2 5(n — 1) T 3ra_4a ’ 


(24) 


dM\ 4t tQ 256t r 3 Q , 2 , 

~dQj s ~~ (ti- 2)T"- 2 _ (3n — 4) Y 3n_4 olJr ^\ a ) 


(25) 


Using Eqs. (MB and m, one can show that the Eqs. (1241) and (l?5l) are equal to Eqs. and (PH), respectively. 
Thus, these quantities satisfy the first law of thermodynamics 


dM = TdS + 


(26) 


A. Stability of the Solutions 


In what follows, we want to investigate the local stability of the charged black hole solutions of Einstein gravity in 
the presence of nonlinear electrodynamics in the canonical and the grand canonical ensembles. In principle the local 
stability can be carried out by finding the determinant of the Hessian matrix of M(Xi) with respect to its extensive 
variables X t , H^ iX = \d 2 M/dXidXj] [24], HH. In our case the mass M is a function of the entropy 5 and the charge 
Q. The number of thermodynamic variables depends on the ensemble that is used. In the canonical ensemble, the 
positivity of the heat capacity Cq = T + / [d 2 M/dS 2 ) is sufficient to ensure the local stability. Since T+ should be 

a positive definite quantity for physical black holes, it is sufficient to check the sign of (d 2 M/dS 2 ) q 


d 2 M 
~dS 2 ) 


\ (n- 2) 


■ (n — 1) i 


+ 


2 (2ti - 3) q 2 


4 (47i — 5) q 4 
■ (n - l) 2 r 5 + n ~ 6 


a + O (a 2 ) . 


Q 


t: (n - l) 2 r^~ 4 


(27) 
















FIG. 8: Asymptotically flat solutions: 10 4 x ^ d g ^l (thin lines) and T+ (bold lines) versus r+ for n = 4, q = 0.01, a 
(solid line), a = 5 x 10 -5 (dotted line) and a = 10 x 10~ 5 (dashed line) 



FIG. 9: Asymptotically flat solutions: ^ ^ (thin lines) and -yj- (bold lines) versus r+ for n = 4, a = 10 4 , and 

(solid line), q = 10.5 (dotted line) and q = 11 (dashed line) 
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Considering Eq. m, we find that the first and second terms are related to the Einstein-Maxwell gravity and 
third one is related to the effect of nonlinearity. In order to find the effects of nonlinearity on the stability of the 
solutions, we plot Figs. [5] and [9] These figures show that there is a lower limit, r+ m i n , for the horizon radius of 
physical black holes (positive temperature). In addition, considering Figs. [5] and [9] one finds large physical black 
holes are not stable. In other words, one can obtain asymptotically flat stable black holes when the horizon radius 
satisfies r+ m i n < r+ < r+ mQX , in which the values of r+ m i n and r+ mQX depend on n, q and a. Although Fig. [8]shows 
that decreasing a leads to increasing r +m i n (slightly increasing r +max ), Fig. [9] indicates that decreasing q leads to 
decreasing both r +min and r +max . 

In the grand canonical ensemble, after some algebraic manipulations, we obtain 


H 


M 

S,Q 



4[(n-l)(n-2)-2 g Vt~ 2TI ] 

(n - 2) r 2 ™-' 2 

1) (n - 2) 2 rX~ 4n - (7 n - 8) q 2 r*+ 6n 
(n — 2) (3 n — 4) 


a 



(28) 


where the last term is the nonlinearity effect of NLED. Regardless of the values of n, q and a, we can write 


ttM I ^ 

S-tQ 1 Small r_|_ 

16 (7 n — 8) q 4 a 
{n - 2) (3 n - 4) r®"" 8 < ’ 

(29) 

jjM I _ 

SiQ 1 Large r. |_ 

4 

(n-l)4"- 2 < ’ 

(30) 


which, in agreement with the canonical ensemble, confirm that the horizon radius of asymptotically flat stable black 
holes should satisfy r+ m i n VP-\-rnax- 


IV. THERMODYNAMICS OF ASYMPTOTICALLY ADS ROTATING BLACK BRANES WITH FLAT 

HORIZON (k = 0) 

Now, we want to endow our spacetime solutions m for k = 0 with global rotation parameters. In order to 
supplement angular momentum to the spacetime, we perform the following rotation boost in the t — <J>i plane 

1 1 —> Et - ai4>i, 4>i i— E(f>i - fet. (31) 

Thus the metric of (n + l)-dimensional asymptotically adS rotating spacetime with p rotation parameters can be 
written as 

/ p \ 2 r 2 p 2 

ds 2 = -f(r) I Edt - ^2 aidfii j + (cadt - l 2 Ed4>i ) + 

\ 2=1 / 2=1 

dr 2 r 2 P o n_1 

— - _ + 7,2 dx ^ ( 32 ) 
^ ' 2=1 2=p+l 

where H = \J\ + 5Zr=i a i/^ 2 - Using Eq. Q, one can show that the suitable gauge potential can be written as 

= h(r) (55° — a,;5^) (no sum on i), (33) 

where h(r) is the same as that in Eq. (JTTJ) . Now, we want to obtain the metric function f(r) for the spacetime (l32l) 
by inserting Eqs. (l32l) and (l33l) into Eq. ©. After some simplifications, we find that the nonzero components of the 
gravitational field equations lead to four different differential equations en, e 22 , 633 and 644 for the unknown function 
/(r), in which 


en = 

ei (static case )| fe=0 , 

(34) 

e 33 = 

e2 (static case )| fe=0 , 

(35) 

e 33 = 

[e2 ( static case) | fc=0 ] (5 2 — 1 ) — l 2 f{r)E 2 [ei ( static case )| fc=0 ], 

(36) 

644 = 

t - f( r ) 1 

v^ 2 l{[ei (static case)\ k=0 ] 9r2 ^ [e 2 (static case )| fc=0 ]}. 

(37) 
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Considering Eqs. m and (EH), we find that the metric function m with k = 0 satisfies all field equations. 
Straightforward calculations confirm that the mentioned rotating spacetime has a curvature singularity at r = 0, 
which may be covered with an event horizon. We can obtain the temperature and the angular velocity of the event 
horizon by analytic continuation of the metric function (1321) and its regularity at the horizon r + . One obtains 


T + (static case ) 


and 

= (39) 

Considering the fact that % = dt + Y^=i ®<t>i * s the nu ll Killing generator of the horizon and using Eq. (fT51) . one 
can find the electric potential as 


, (38) 

k= 


^ (static case) 


(40) 


Here, we calculate other conserved and thermodynamic quantities of the black brane solutions. Like previous section 
and with the same approaches, one can show that the entropy and the electric charge per unit volume V n -\ of the 
presented black branes are, respectively, given by 


and 


S = 


ErT 1 



(41) 


(42) 


Now, we should calculate the finite mass. In general, the action Iq, diverges when evaluated on the solutions, as 
the Hamiltonian and other associated conserved quantities. To compute the conserved charges of the asymptotically 
adS solutions of Einstein gravity, we use the counterterm method S3- This method was inspired by the anti-de 
Sitter/conformal field theory (AdS/CFT) correspondence and consists in adding suitable counterterm I ct to the 
action Iq in order to ensure the finiteness of the boundary stress tensor derived by the quasilocal energy definition 
[28j. For asymptotically adS solutions of Einstein gravity with flat boundary, the suitable counterterm I ct is given by 


ht ~ 8 7T , 


d n Xyf^ 


'dM 


n — 1 
l 


(43) 


Varying the total action (Itot = Ig + let ) with respect to the induced metric 7 ab , we find the boundary stress-tensor 


as 


rpab _ 


87T 


0 Qb - (0 


n — 1 


(44) 


Now, we choose a spacelike surface B in <9A4 with metric ay, , and write the boundary metric in Arnowitt-Deser- 
Misner form 


'y lll/ dx ,1 dx v = —N 2 dt 2 + <Tij (dip 1 + V l dt) (dp? + V^dt) , (45) 

where the coordinates p l are the angular variables parameterizing the hypersurface of constant r around the origin, 
and N and V" 1 are the lapse and shift functions, respectively. When there is a Killing vector field £ on the boundary, 
then the quasilocal conserved quantities associated with the stress energy momentum tensor of Eq. (PI) can be 
calculated as 


Q (0 = f d n ~ 2 p^T ab n a t; b , (46) 

JB 

where a is the determinant of the metric ay,-, and n a is the timelike unit normal vector to the boundary B. The 
conserved quantities associated to the timelike £ = d t and rotational ( = d^ Killing vector fields are 

M = J d n ~ 2 p^T ab n a i b = (nE 2 - l) , (47) 
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FIG. 10: Asymptotically adS solutions: y 9 d J$ J (thin lines) and -y±- (bold lines) versus r+ for n = 4, q = 100, S = 1.1 and 
A = —1, and a = 0.001 (solid line), a = 0.02 (dotted line) and a = 0.05 (dashed line) 


Ji = [ d n 2 py/aT ab n a ( b = 

Jb 16 tt 


(48) 


which are the mass and the angular momentum of the system enclosed by the boundary B. To check the first law of 
thermodynamics, we obtain the total mass M as a Smarr-type formula 


M (S, J, Q) = 


(■nZ - 1) J 


rds/Z{Z- 1 )’ 

where Ji = JT J 2 and Z = S 2 is the positive real root of the following equation 

2AT- 2 167TJ4'"- 2 32t r 2 Q 2 4' 2Tl " 4 1024 7 r 4 Q 4 4' 4n - 6 


• (n — 1) nlyJZ (Z - 1) (n—l)(n — 2)Z [n - 1) (3n - 4) Z 2 


a + O (a 2 ) = 0, 


(49) 


(50) 


with 4/ = 


Vz{4,sy 


l/(n—1) 


. It is a matter of straightforward calculation to show that the conserved and thermo¬ 
dynamic quantities satisfy the first law of thermodynamics 

dM = TdS + Y VidJ z + (51) 

i 

In other words, the quantities T = (5^) and $ = are the same as those calculated in 

Eqs. (f38|) , (1391) and (l40l) . respectively. 


A. Stability of the Solutions 


The final step is devoted to analyzing the local stability of charged rotating black brane solutions of Einstein gravity 
in the presence of nonlinear electrodynamics. We use the similar theoretical manner as was discussed in the Sec. m 
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FIG. 11: Asymptotically adS solutions: ^ (thin lines) and Sj- (bold lines) versus r. f. for n = 4, q = 100 and A = —1, 

a = ICG 5 with H = 1.1 (solid line), S = 1.2 (dotted line) and 3 = 1.3 (dashed line) 


and investigate thermal stability in both the canonical and the grand canonical ensembles. In the canonical ensemble, 
the electric charge and the angular momenta are fixed parameters, and (d 2 M/dS 2 ) 7 q at constant charge and angular 
momenta is 


fd 2 M\ _Ci C 3 
V dS 2 ) jq C 2 + C 4 a + ° (a) 


Ci = —2 {n [3 (n — 2 )^ 2 — n + 3] q 4 r + 3n+6 + 
2 [3(n - 2)S 2 - n 2 + 3] Ag 2 r +"+ 4 + 

(n - 2 ) [(n + 2 ) 2 2 - n - l] A 2 r™+ 2 } , 


C 2 = (n- 1 ) 7 rS 2 [(n - 2)2 2 + l] [(n - 2)Ar+" - nq 2 r \] , 


C 3 = -4g 4 {n 2 [(17n 2 - 54n + 40)S 2 - n(5n - 23) - 20] g 4 r+ 5n+10 - 
2n(n - 2) [(13n 2 - 50n + 40)S 2 - n(n - 19) - 20] A q 2 r+ 3n+8 + 
[n - 2) 2 [(5n 2 - 42n + 40)S 2 + n(7n + 11) - 20] A 2 r+ n+6 } , 


(52) 


(53) 

(54) 


(55) 


C 4 = (n — 1)(3 n — 4)n^ 2 [(n — 2)^, 2 + l] [(n — 2)Ar+" — nq 2 r+] 2 , (56) 

where in the Eq. (1521) the first term is related to the Einstein-Maxwell gravity. 

Here, we plot Figs. [TO] and [FT] to investigate the nonlinearity as well as rotation effects. These figures show that 
for suitable fixed values of the metric parameters, there is an r +m j n , in which for r+ > r+ m j n we can obtain physical 
asymptotically adS rotating stable black brane solutions. In other words, we find that, unlike the asymptotically flat 
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solutions, the asymptotically adS rotating large black brane solutions are stable. Considering Fig. 1111 one may find 
the rotation parameter affects on the values of both T + and 

In the grand canonical ensemble, one can show that the determinant of the Hessian matrix has the following form 


ttM 

U S,Q,J 


647T ( q 2 - A 2 ) 

B[ 


Bo 


Bi(3n — 4) [ nq 2 — (n — 2)Ar 2n 2 ] r 


+ O (a 2 ) 


(57) 


Hi = l 2 Z 6 [(n - 2) S 2 + 1] [nq 2 - (n - 2)A r 2 + n ~ 2 ] 


2n— 2"| 3n—4 

r + ’ 


(58) 


1?2 = 2567t< 7 2 [3(n — 2) 2 A 2 r+" 4 + 3n(n — l)q 4 — 2(n — 2)(3n — 2)Aq 2 r 2 [ n 2 ] , (59) 

where in the Eq. (1571) the first term is the determinant of Hessian matrix of Einstein-Maxwell gravity and the last 
term indicates the nonlinearity effect. Following the method of previous section and regardless of the values of n, q , 
A, S and a, one finds 


IjM I 

S,Q,J I Small r_|_ 


I 

S,Q,J I Large r_|_ 


7687r(n — 1 )q 2 a 


n(3n — 4)Z 2 ^ 6 [(n — 2) S 2 . _ 
647T 

(n - 2)Z 2 S 6 [(n - 2) S 2 + 1] r^" 4 


_5n—6 


<o, 


>o, 


(60) 

(61) 


which are in agreement with results of the canonical ensemble and confirm that large black branes are stable. 

In order to analyze the correctness of our discussions for stability criterion of asymptotically flat and adS black 
objects, we should argue for the validity of numerical calculations. Regarding Eq. m, we find that higher order 

/ \ ( 21 + 1 ) 

terms of electric field can be formed by increasing j in I 9 _ 1 J (j = 0,1, 2,...). Therefore, ignoring the higher 

( \ {2j+1) 

order terms of a makes sense, if increasing j leads to reasonable decreasing of ( I ad . In the following 

mi—i ^ r+ ' 

tables, we consider the numerical calculations of Figs. 181 fill 


Fig. E 

Fig. [9] Fig. [10] 

Fig. HD 

q = 0.010 

q = 10.00 q = 100 

q = 100 

a = 10 -5 

a = 10 -5 a = 10 -3 

a = 10- 5 

r + « 0.075 

r+ ss 2.40 r+ ~ 4.6 

r+ S3 4.64 

4 ~ 23.704 

r + 

4 ~ 0.72 4 ~ 1.03 

r + r + 

* - 100 

0.133 

a S3 0.38 x 10 -5 a S3 10~ 3 

(4a«10- 3 

(yr)° a 2 ~ °- 75 x 10 ~ 3 

(^-)° a 2 S3 0.20 x 10 _1 ° (^)°0! 2 S3 !°- 6 

(t)V»io-‘“ 


Left table: asymptotically flat solutions & Right table: asymptotically adS solutions. 


Taking into account the numerical results of the tables, one can confirm that numerical calculations of stability 
conditions for both asymptotically flat and adS black objects are logical. 


V. CONCLUSIONS 

Motivated by the (quartic) string corrections of Maxwell field strength, at first, we obtained black hole solutions 
of Einstein-NLED gravity with various horizon topology and investigated their geometric properties. Then, we fixed 
A = 0 and k = 1 to calculate the conserved quantities of the asymptotically flat black holes. We obtained a Smarr-type 
formula for the mass as a function of the entropy and the electric charge of the solutions and checked the first law of 
thermodynamics. We studied the stability analysis of the asymptotically flat black holes both in the canonical and 
the grand canonical ensembles and investigated the effects of NLED. We found that for the fixed values of n, q and 
a, small and large physical black holes are not stable. It means that obtained asymptotically flat black holes can be 
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stable when the horizon radius satisfies r+min < r+ < r+ max , in which the values of r+ m in and r+ m ax depend on 
the metric parameters. In addition, we found that although decreasing the nonlinearity parameter leads to increasing 
r+min (slightly increasing r+ m ax)i decreasing the charge parameter leads to decreasing both r +m i n and r+ m ax- 

After that, we considered the horizon-flat solutions and used a suitable rotation boost to endow angular momentum 
to the asymptotically adS spacetime. Using the counterterm method, we obtained the conserved quantities of the 
asymptotically adS black branes. We also obtained a Smarr-type formula for the finite mass as a function of the other 
quantities and showed that they satisfy the first law of thermodynamics. Besides, we performed a stability analysis 
of the rotating solutions both in the canonical and the grand canonical ensembles. We showed that there is a lower 
limit, r+ m i n , for the physical solutions (positive temperature). Stability analysis of both ensembles confirmed that, 
unlike the asymptotically flat solutions with spherical horizon, the horizon-flat asymptotically adS rotating black 
brane solutions with large event horizon are stable. In other words, we showed that there is an r+ min for suitable 
fixed values of the metric parameters, in which for r+ > r+ m i n , the asymptotically adS rotating black brane solutions 
are stable. Moreover, we fixed the values of n, A, a and q to analyze the rotation’s effect on the stability conditions. 
We showed that although 2 does not change the location of r +m i n , it can change the values of the temperature, the 
heat capacity and the determinant of Hessian matrix. 

It is notable that due to the negative temperature, the small black holes/branes are not physical and we should 
restrict the horizon radius to r+ > r+ m in, while (in)stability of large ones is related to their horizon geometries. In 
other words, large black holes (branes) with k = 1 (k = 0) are unstable (stable). 

Finally, it was seen that the nonlinearity part not only modified electromagnetic part of the solutions, but also the 
kind of horizons and thermodynamics properties. In absence of of correction part, with suitable choices of parameters, 
metric function could acquire two horizons whereas, for this nonlinear theory, it is possible to find three horizons. This 
fact has some application regarding anti-evaporation of black holes/branes. The structure of black hole in presence 
of nonlinear electrodynamics is quite different comparing to the linear Maxwell theory and its phenomenology is 
also describing a more general case. In addition, it is worthwhile to mention that considering the first order effects 
of nonlinear electrodynamics changed the properties of the black objects at small distances. In other words, this 
generalization changed timelike singularity of Reissner-Nordstrom black holes to spacelike singularity. Hence, in order 
to recover the properties of the Reissner-Nordstrom black holes, it may be logical to keep terms only upto quadratic 
order of a in the series expansions. Besides, it is notable that one can consider asymptotically adS black holes with 
spherical (k = 1) and hyperbolic topologies [k = —1), to investigate P — V criticality in the extended phase space of 
the solutions by calculating the Gibbs free energy for various A. These extensions are under examination. 
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